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Abstract. We present a theoretical study of superconductivity of polarons in the Hubbard-Holstein model. 
A residual kinematic interaction proportional to the square of the polaron hopping energy between polarons 
and phonons provides a pairing field for the polarons. We find that superconducting instability in the d- 
wave channel is possible with small transition temperatures which is maximum in the large to small polaron 
cross-over region. An s-wave instability is found to be not possible when the effective on-site interaction 
between polarons is repulsive. 

PACS. 63.20.Kr Phonon-electron and Phonon-Phonon interactions - 74.20.Mn Polarons and bipolarons 
- 74.20.Fg BCS theoy and its development 

1 Introduction interacting with local phonons is the Hubbard-Holstein 

model [lU^. If the e-ph interaction is large enough, it can 
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Polaronic superconductivity has been a subject of interest jead to the formation of polarons. Building on an earlier 

• I— I 

|Xj following discovery of superconductivity in cuprate oxides ^o^j, ^ ^ Xakada and Hotta recently studied [5] this model 

. . . . and some molecular conductors Ij. A characteristic fea- ^j^g strong e-ph couphng hmit where the polarons are 

ture of these materials is that the conduction electrons gmall polarons. For intermediate e-ph coupling, however, 

live in a narrow energy band. It is also recognized that ^he polarons may not be small polarons depending on e- 

electron-phonon (e-ph) interaction is considerable in these coupling strength and the ratio of the values of the 

materials. The simplest model for narrow band electrons electronic hopping to the phonon energy. In this paper we 
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the role of large-to-small polaron crossover on supercon- through its maximum in the large to small polaron cross- 
ductivity. over region (see Fig. 1). The details of our calculations 

are given in the rest of this section and the conclusions 
are given in Sec. III. 

The application of MLF transformation to H leads to 
Hmlf = e^He~^, where R is given by, 



2 Mean-field theory of superconductivity in 
the Hubbard-Holstein model 



The Hubbard-Holstein (HH) model [H is: i? = Ao ^ ~ + Ai ^ - b,+s) , (2) 

H = — t cj^c^-|^(5a- + [/^^nj|ni|^ Aq and Ai represent the lattice deformations at the elec- 

iSa i 

+ J2^^b^+9Y.'^Ubl+b^)■ (1) 



6 Ci(jC — Ci(j CXp 



5 

(3) 

and, 



tron site and its nearest-neighbor sites, respectively. They 

* are treated as variational parameters to be determined by 

Here t (> 0) is the hopping energy between molecules at ..... / -.i ^ ^ \ J\^r j^j. 

^ ' ± to bj mmimization (with respect to Ao and Ai) of ground state 

lattice site i and its nearest-neighbor lattice sites i -\- 5, r ii j- r j tt -li • tt • ii u 

' energy oi the transiormed Hamiltonian. Using the above 

Cicr (cjg.) is the annihilation (creation) operator for the obtain' 
electron with spin a at a lattice site i and riicr is the cor- 
responding number operator, U is on-site Coulomb repul- 
sion, bi (hi) is the phonon annihilation (creation) oper- 
ator, and g is the e-ph interaction strength. All energies 

are measured in terms of the phonon energy {hui). In the ^ — bi — Xo riia- — Ai rii+s.cr ■ (4) 

calculations given below, first we apply a Modified-Lang- r^^^ transformed Hamiltonian (Hmlf) is obtained as: 
Firsov (MLF) transformation 6 (which is the original LF 

Hmlf = Ho + Hi + H2 , (5) 

[7] transformation modified to include lattice deformations 

on the sites nearest neighbor to the site where an electron where, 

resides) to H to convert it to the polaron representation, tt _ ST^ uiu _ , „ , rr „ 

Ha — ^ ^ ^i^i 2_. ia ~r ^ ef f / ^ 

then we eliminate the residual polaron-phonon (pol-ph) * 

interactions to obtain an effective interaction between po- «<5 <^<^' 

+ V2 y ni+s,crni+s' ■, (6) 

larons, and finally we do a mean- field theory of supercon- i5s^=£S' aa' 

ductivity specializing to the cases of s and d wave order = -t ^ c|'^Cj+5,<^ exp{X^ ~ X.+a) , (7) 

parameters in the case of a square lattice. Our main result 

and 

is contained in Fig. 2, where it is shown that the supercon- 

A . . • fi. ^ 1. 1 = (ff-Ao)^n,„(6l + &.) 

ducting transition temperature m the a- wave channel goes 
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- Ai ^ nia{bl_^_s + bi+5) . (8) where, 

Here, H3 = -tp ^ cl^Ci+s,a [(Ao - Ai)(6t - 6^ - bl^^ + bi+s) 

€p = Xo{2g -X,)-zXi, (9) + Ai ^ (6]+^, - bi+s') 

U^ff = U-2ep, (10) - ^1 (^l+s+S' - bi+s+5')] (14) 

Vi =2(5-Ao)Ai,y2 = A?, (11) 



5' 5^-5 

and, 



Xi = Ao(6t - bi) + Ai J2ibl+s - bi+s) , (12) 



tp = t exp [-(Ao -Xif-{z- 1)A?] . (15) 



5 

and z is the coordination number of a lattice site. The 



H = Ho-tpY^ clc^+s,a + Hi, (16) 



i6cr 



In the above equation H2 and describe the pol-ph inter- 
actions. On eliminating the pol-ph interaction terms and 

Eq. (5) is the Hamiltonian of polarons. The polaron is . j j 1 uj. • /r 

^ ^ ' ^ ^ considering second order processes only, we obtain (for 

much more massive than the original electron since it , ^ 

ztp < 1), 

is a composite particle consisting of an electron and the 
associated lattice deformations. These polarons interact 
among themselves (through the interactions Ugff, Vi and -^^iiere 
V2), and interact with the phonons through off-diagonal 

H4 = -2{g - Ao)^ ni-ffiii 

phonon terms of Hi and H2. The phonons remain unaf- 

fected by the MLF transformation. Notice also that the + 2A,(,-Ao)EE 

iS <t<t' 

pol-ph couplings are much weaker than the original e-ph — ^ Ui^s <T'^i-\-6' a' 

iSS' aa' 

coupling (51). 2/\ t 

i55' aa' 



Now we work using the variational phonon basis ob- "^"^ 



tained through the MLF transformation. The variational 
parameters Ao and Ai are found out from a minimization 



- 4a'4a(^i+S,'7Ci+S',a']- (17) 



of the ground-state energy in the normal state with zero It may be noted that our study is restricted to the region 

phonon averaging. The values of Aq and Ai depend on t/uj where tp{Xo - Ai), (g - Ao), and Ai are much less than 

and g/uj for a particular lattice and not on U. With this 1 even though g > 1. This justifies our perturbation cal- 

variational basis we proceed with our second order pertur- culation. On simphfication, confining to nearest neighbor 

bation calculation by expanding the exponentials in Hi to singlet pairing terms, and neglecting terms 0(A?) since Ai 

obtain, up to 0(Ao - Ai), H valid for (Ao - Ai) < 1 as: is quite small (see Fig. 2), we get: 

H = Ho - tp^cl^Ci+5^a + H2 + Hs, (13) H = -tp^cl^a+s^a - ep^riia + U^ni^nii 

i6<T iS<T i<7 i 
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tliXl - 2AoAi) XmK>5,^'4^c,+5,^c„ 



=< 4+sAl > ' (21) 



(18) 



where U = {U — 2g^). The elimination of the pol-ph terms 
generates a pairing field whose magnitude is 0{t^) which 
is rather small. Note that the field associated with the on- 
site pair-hopping processes (the last two terms in H) is 
repulsive. The origin of superconductivity in this model, 
when [/ > 0, then is the field generated through a virtual 
emission and absorption of phonons while the polarons 
hop from site to site. The two-polaron process involved is 
one polaron hops to its nearest neighbor site and excites 
a phonon at the target site and another polaron on the 
target site absorbs the excited phonon and hops to the 
site left by the first polaron. These processes lead to the 
fourth and fifth terms in If the e-ph interaction is not suf- 
ficient enough to lead to the formation of polarons, super- 
conductivity occurs in the usual way through the phonon 
exchange between electrons, and the pairing field contains 
the full e-ph interaction. Now we are at a stage where we 
can do a mean-field theory of the superconductivity using 
H. 

On mean-field factorization of H , we obtain: 

Hmf = -tp ^ -i'^p-Ui^)Y^ " 

+ [/^(AoC^CT + h.c.) - iVc7|Ao|^ 

i 

+ Vtp Y^{Asc^+sic,^ + h.c.) - NzVtp\As\^ , (19) 



where 



and 



U = U-zVtp, 



Vtp = -U'p{Xi - 2AoAi) 



(22) 



(23) 



The above mean-field approximation is valid for U less 
than the polaron band-width. For larger U strong corre- 
lation effects like Hubbard sub-band formation can occur. 
For recent studies of strong correlation effects originating 
from large U in Hubbard-Holstein model see Refs. [HIH]- 
Converting to momentum space, we get: 



Un 



E 



kcr \ / kcr 

UAo + VtpAi{k) ] c_kiCkT + h.c. 



where 



NzVtp\Ai\^ 



C(k)=^e*^ 



Ai(k)=^zi 



56 



ik.5 



(24) 



(25) 



(26) 



(20) 



and N is the number of lattice sites. In the s-wave case, 
U term should be considered even if J7 > since the ex- 
tended s-wave order parameter has an on-site component, 
and this is suppressed by the on-site repulsion V . The 
V term does not affect c?-wave pairing. In studying su- 
perconductivity, we will be considering only ?7 > case 
since J7 < case leads to an effective attractive Hubbard 
model situation and this model is well understood pUlfTTl 
112] . Now we will consider s and d wave pairing channels 
separately. 
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2.1 s-wave channel 



In the s-wave case, the effective Hamiltonian is 



as 



-ff s = ^ e(k)4^Ck<T + ^(Ds(k)c_kiCkT + h.c.) + So , 

ka k 

(27) 

where, ?(k) = -tp^{k) - /i, D,(k) = UA, + FtpAf(k), 
Af(k) = 2Z\57s(k), = n + ep-Un/2, 7,(fc) = cos(fc:ra) + 
cos{kya) for a square lattice of lattice constant a, and 
s„ = -(iV[7n74) - nU\A^\^ - NzVtp\As\^. In the above 
we have introduced a chemical potential (/i) to fix the 
number density, and ^i(fc) is for the isotropic case. The 
integral equation which determines I?s(k) is obtained, us- 
ing Green's functions, to be. 



DJk) = - — Vv^JkkO ^'^^'\ tanh 
^ ' N ^ '2Es{k') 



where the pairing field ys(kk') is 



F,(kk') = U + 4z-Wtplsikhsik'), 



(28) 



(29) 



/3 — l/fc^T, ks is the Boltzmann constant, and the quasi- 
particle energy Eg (k) = y^e'^{k) + \Ds{k)\'^. Now, we have 
already mentioned that we will consider the non-trivial 
case oi U > 0. In that case s-wave pairing is not possible 
since T4(kk') is not attractive. 

2.2 d-wave channel 

In this case, we have 

-kiCkT + h.c.) + do , 



k(T 



(30) 



PEd{k') 



(31) 



where, 

Vd{kk') = 4z- Vtp7rf(k)7rf(k') , (32) 

Ej{k) = e2(k) + |L'd(k)|2, and7<j(A:) = cos{kxa)-cos{kya) 
for a 2D-square-lattice case. The superconducting transi- 
tion temperature is obtained from. 



k' 



2e~(k') 



tanh 



(33) 



where £'d(k) 



VtpAi{k), and do 



-{NUn^/A) 



where /3c = l/ksT^f 

We have numerically analyzed the above equations and 
the results are shown in Figs. 1-2. The values of t (< 0.5) 
used are in the anti-adiabatic limit. In Fig. 1 we have 
displayed the variation of the lattice deformation (Ao and 
Ai) as a function of g. A smooth cross-over from a large 
to a small polaron 6J with increasing g is seen in Fig. 1. 

The variation of Tc with g (for n — 0.81) for supercon- 
ducting instability in the d-wave channel is displayed in 
Fig. 2. The variation of Tc with n has a bell shaped form 
centered at half-filling similar to that found in Rcf. [S]. 
Superconducting instability with small is possible for 
intermediate values of e-ph coupling. The small is due 
to small value of the pairing field [which is 0{tp)]. The 
is seen to reach a maximum in the large to small polaron 
cross-over region. 

3 Conclusions 

In this work we presented a study of superconductivity 



NzVtp\As\'^. Then the c?-wave gap equation is obtained in the anti-adiabatic limit of the Hubbard-Holstein model 
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considering the polaron formation. Our purpose was to 
study superconductivity of the polarons in the large to 
small polaron cross over region. We found that transition 
temperature for superconducting instability in the d-wave 
channel goes through its maximum in the large to small 
polaron cross-over regime. We also found that s-wave pair- 
ing is not possible for the case of U > 2g^. In the strong 
e-ph coupling limit also, it was found [5j that is ex- 
tremely small when U > 2g^. In real materials there is 
inter-site Coulomb repulsion and this will drastically re- 
duce [13] the T^. So, we conclude that, when the mech- 
anism considered here operates, is very small both in 
the intermediate and strong coupling limits of the e-ph 
interaction. 
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Fig. 2. vs g for n = 0.81 and t: 0.35 (thick line) and 0.30 
(thin line). 



